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■ Two-body non-lcptonic charmed decays £>( s ) — > D^P, D^P, D( S )V an d D^^V are 

analyzed in perturbative QCD approach, where P and V denote the light pseudoscalar 
meson and vector meson, respectively. We test the D meson wave function by a x 2 fit with 



(N 



experimental data of six B — > DP channels. We give the branching ratios of all the charmed 
B decay channels, most of which agree with experiments amazingly well. The predicted B s 
decays can be confronted with the future experimental data. By straightforward calculations, 
' our pQCD approach gives the right relative strong phase of 02/01 with experiments. We 

^ also predict the percentage of transverse polarizations in B^ — > D*V decay channels. 
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I. INTRODUCTION 

B physics experiments provide a good test of the standard model and severe constraints of 
new physics parameters. Recent years more and more efforts have been made to the study of B 



meson decays both experimentally and theoretically. In the near future, there will be more and 

o ; 

more data in B physics thanks to the run of B factories, Tevatron and LHCb. Theoretically, a 



great improvement 



ras been made to the study of exclusive decays of B mesons. In the history, 



s Q], t 



Currently, perturbative QCD 



^ . naive factorization 

failed to explain color-suppressedprocesses such as B° — > D° 
factorization approach (pQCD) \$ is one of the popular methods to deal with the two-body non- 
leptonic decays of B mesons. It explains the experiments successfully, especially for the direct CP 
asymmetries when the final states are two light mesons, which inspires people to see how far it 
will go. 

Charmed decays of B {B s ) mesons are more complicated than the decays with only light mesons 
as final states. The B — > D transitions involve three scales: Mb, Mp, A, respectively. The 



factorization was proved in soft-collinear effective theory [6j with less predictions than pQCD 
approach, since they need more inputs than pQCD approach. In B — > light transition, the light 
spectator quark in B meson is soft, while it is collinear in the final state meson, so that a hard gluon 



2 




FIG. 1: color allowed diagrams in pQCD approach for B — ► DP decays 



is needed to connect it to the four quark operator shown in Figure 1. The momentum square of the 
hard gluon connecting the spectator quark in B — > D transition is only a factor of (1 — m 2 D /m B ) 
than the B — > light transition, which shows pQCD should also work well in B — > D transitions. 
The hierarchy: Mg » Mp >> A is used in pQCD framework [jj]. Some separate calculations on 
B — > D decays in pQCD approach are carried on [§[ in the leading order of Mb /Mb and A/ Mb 

expansions. It's found that the pQCD do work well since the D meson recoils fast. 

(*) 

In this paper, we calculate all the processes of a B^ meson decays to a D^j meson and a light 
pseudoscalar meson or vector meson. Only tree diagrams contribute to these processes involving 
only one kind of CKM matrix elements. So there's no direct CP asymmetry in these decays. The 
light cone distribution amplitudes (LCDA) of mesons are necessary inputs in the pQCD framework. 
The light mesons' distribution amplitudes (DAs) have been well studied and calculated by QCD 
sum rules. But there are little studies on heavy mesons' DAs, especially for D meson. In this 
paper we collect several candidate distribution amplitude models for D meson, then we fit out the 
parameters using the experimental results and make a comparison among them. 

The paper is organized as follows: Section [II] contains the conventions and notations that we 
adopt, together with all the wave functions used in this paper. The pQCD analytic formulae for 
the amplitudes are given in section IIIII Section IIVI contain the numerical results and discussions. 
Section [V] is our summary. 



II. ANALYTIC FORMULAE 



For the charmed B decays we considered, only the tree operators of the standard effective weak 
Hamiltonian contribute. The Hamiltonian is given by: 

n eff = %v cb v: q , [Cxi^oM + c 2 ( m )o 2 (m)] . (i) 
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The tree operators are 

Oi = {c a bp)v-A{q' pU a ) V - A ,0 2 = (c a b a ) V _ A (q'i3Uf3)v-A , (2) 

with q' = d,s. a (3 are the color indexes, and (qiq 2 )v-A = 9i7^(l — 7 5 )<Z2- The V qiq2 in the 
Hamiltonian denote the CKM matrix elements. 

In the factorization formulae, the combinations of Wilson coefficients usually appear. In this 
paper, they are defined as follows: 

a x = C 2 + Ci/3 , a 2 = Ci + C 2 /3 . (3) 

In this paper, all the momentum are defined with the light cone coordinate. Two light-like 
vectors n and v are defined, with n 2 = 0, v 2 = and n ■ v = 1. n can be written as (1, 0, Ot), and 
v is (0, 1, Ot)- The momentum of B meson, D meson and the lightest meson are denoted as Pi, P 2 
and P3 respectively. At the rest frame of B meson, the light meson moves very fast. So or P^ 
can be treated as zero. In this paper, the momentum are defined as: 

Pl = (1,1, Or) , P2 = (1^0 T ) , P 3 = ^(0,1 - r'.Or) . (4) 

The momentum of the light anti-quark in the B meson and D meson are denoted as k\ and k 2 
respectively, while ks is defined as the quark momentum of lightest meson. In our calculation they 
are taken as: 

u l Mb n u ^ h /n (1 ~ r 2 )M B 

= 1/2 ' 3 = ^ ,X3 7f ' k 3±j, 

M B 

ki = (x\—=-, 0, ku_) for color suppressed contributions, 
V2 

M B 

k\ = (0, xi—=, ku_) for the others. (5) 
y2 

with x\, x 2 and X3 as the momentum fraction, and ku_, k 2 j_ and k3j_ as the transverse momentum 
of the quark. 



A. Wave functions of Bf s \ mesons 

In pQCD calculation, the light-cone wave functions of the mesons are needed. The B meson 
and B s meson have the similar structure of wave function, except different values of parameters 
characterizing a small SU(3) breaking effect. In general, the Br s ) meson wave function are always 



4 



decomposed into the following Lorentz structures: 

r e*-*(0M0Hj(*)|B(.)(i\)> 



(2vr)4 

There are two distribution amplitudes in the above equation. However, 4>B^(ki) gives smaller 
contribution We will neglect it in our calculation and only keep the first Lorentz structure. 

*b w = ^=(/Pi+M B(3) ) 75 <j> Bis) (k 1 ). (7) 



The distribution amplitude in the b-space is: 



<f>B (s) (x,b) = N B( s) x 2 (l - x) 2 exp 



(8) 



2 wj 2 

with b as the conjugate space coordinate of ku_. Nb, s s is the normalization constant, which is 
determined by the normalization condition: 

dx( f> B(s) (x,b = 0) = ^=. (9) 

For parameter ujf,, we take the value 0.40 GeV for and mesons, and 0.50 ± 0.05 GeV for B® 



meson, characterizing the small SU(3) breaking effect as argued in 



B. Wave functions of light pseudoscalar mesons 

The decay constant of the pseudoscalar mesons is defined by: 

<0|gi7„75<?2|P(P 3 )} = ifpP». (10) 

The decay constant for ir and K are /„• = 131 MeV, fx = 160 MeV. 

The light cone distribution amplitudes (for out-going state) for light pseudoscalar mesons is: 

(P(P 3 )\q 2(3 (z) qia (0)\0) (11) 

JxP 



dxe™ p - z [75 (x) +75m o (x) + m j 5 (j6 jh - l)4> (x)] , 

where u is the light cone direction that the light pseudoscalar meson's momentum is defined on. 

M 2 

The chiral scale parameter tuq is defined as uiq — — 



m 0l +m„ 



5 



The distribution amplitudes are usually expanded by the Gegenbauer polynomials. Their ex- 
pressions are: 



3/p 
fp 



x(l - x) 1 + afcf 2 (i) + a$C?\t) + a£(7f (t) 



^3/2, 



2 v / 27Vc 
/p 



l + a p 2 C l 2 l2 {t) + alCl ,2 {t) 



with t=2x-l. The coefficients of the Gegenbauer polynomials are 



2® 



at = 0.44 , at = 0.25 , af K = 0.17 , a$ K = 0.2 , 
a p 2w = 0.43 , = 0.09 , a\ K = 0.24 , a\ K = -0.11 



a 3lT = 0.55 , a 



3K 



0.35 



(12) 
(13) 
(14) 



(15) 



For 77 and 7/, the mixing mechanism must be taken into consideration. We will take the method 
presented in ref. [l2|. where the r] n and rj s are chosen as the basis of mixing: 



U(</>) 



(16) 



with 



—={uu + dd) , |r/ s ) = ss , 
V2 



(17) 



C/(0) 



COS 1 

sin< 



- smi 
cos c 



(18) 



where the mixing angle cj> = 39.3° ± 1.0°. 

The assumption that the distribution amplitudes of r\ n and r\ s is the same as the distribution 
amplitudes of ir is adopted, except different decay constants and chiral parameters. The decay 
constant of r] n and r/ s is taken from ref. [3]: 

f n = (1.07 ± 0.02)/; = (139.1 ± 2.6) MeV, f s = (1.34 ± 0.06)/^ = (174.2 ± 7.8) MeV. (19) 



And the chiral parameters are given by 

1 



TOn 



- \rni cos 2 6 + m 2 ., sin 2 c 

r 2 2 1 , 2-2 

[m 1 cos <p + to sm 1 

2to s ' ' 



v2/s / 2 2 \ j. • 

— — (m^/ — m^) cos sm q 

In 

{l! 1 (to 2 / — to 2 ) cos 6 sin < 
V2fs V n 



(20) 
(21) 



TABLE I: The decay constants of vector mesons (in McV) 

fp flC* fui f<j> fj S'k* fu ftp 

209 ± 2 217 ± 5 195 ± 3 231 ± 4 165 ± 9 185 ± 10 151 ±9 186 ± 9 



C. Wave functions of light vector mesons 

The decay constants for the vector mesons are defined by: 

(0\qilv<l2 1 V(p, e)) = fvmve^, (O^a^lVip, e)) = ify {e^p u - e u p^). (22) 

The longitudinal decay constant of vector mesons can be extracted experimentally [jij. And the 
transverse ones can be calculated by the QCD sum rule [bJ ]. We list all of them in Table HI 
The distribution amplitudes up to twist-3 of vector mesons are 



(V(P, el)\q 3P (z)q la (0)\0) 



l = J dxe lxP - z [M V AUv(x)+ Al + Mv4> s v(x)] a/3 , 



(V(P,e* T )\q 2p (z)q la (0)\0) = J dxe lxP - z [M v ji^ v v (x)+ 

+M v ie^ pal ^e* T »nPv a ci ) a v {x)} a(] , 



(23) 



where x is the momentum fraction of the q 2 quark. And convention e 0123 = 1 is adopted for the 
Levi-Civita tensor. The twist-2 distribution amplitudes of vector mesons are expanded as: 

3/y 



<f>v(x) 



3/v 



x(l — x) 



x(l — x) 



%/2iVc 

We take the following values for the Gegenbauer moments 



l + a\cl(t) + 4cl(t) 



l + aj;C?(t) + a 2 L C 2 2 (t) 



(24) 



25|: 



0.15 ± 0.07 , a\ K * = 0.03 ± 0.02 , a% K . = 0.11 ± 0.09 , = 0.18 ± 0.08 



a£ = = 0.14 ± 0.06 , a^ K * = 0.04 ± 0.03 , a£ K . = 0.10 ± 0.08 , aL = 0.14 ± 0.07 . (25) 



For the other distribution amplitudes, we use the asymptotic form: 



3f v (-t) 



3/y 



4^6 



(26) 
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D. Wave function of £)W meson 

The two-particle light-cone distribution amplitudes of meson, up to twist-3 accuracy, are 
defined by [7|: 

'■™(0\cp(0)d 7 (w)\D*-(P)) = -—L^[(f + M D .) M L D ,{x) + {P + M D ,) ji T 4>l*(x)U 



I 



(2vr) 4 / n y| v " y/2N t 



c 



(27) 



with 



dx<f> D (x) = — ^= , [ dx<f>k*(x) = , [ dx(j) T D *{x) = -lm= (28) 

as the normalization conditions. In the heavy quark limit we have 

T m c + m d rp m c + m d n( T lM s (0Q s 

Id* ~ Jd* M ~ Jd* ~ Jd* m ~ U(A/Md*). (29) 

Thus we use = fo* in the calculation. There are several candidate distribution amplitude 
models for D meson. We collect as below: 

A Gen) M) = ^=fD6x(l-x)[l + C D (l-2x)], 

2 h 2 



1 . „ ^ . „ r -UJ Z b 



^ Gen) (x,b) = _^= /D 6x(l-x)[l + C D (l-2x)]exp[ 

^ = ^=^iV D ^(l - x) ex P [- - (— -) -_], 

-(GOV), ,x 1 , AT r xMp 1 

n (£, 0) = — ^= f niVnx exp — k, 



^ uan9 \x,b) 



1 „ , T \ r . (1 — x)m~, + xm 2 , . 

— /oJVDx(1 - x )exp[-A ^ (30) 



In the above models x is the momentum fraction of the light quark in D meson. The first DA 
model (p^ en ^ was proposed in Q], which is the Gegenbauer polynomial-like form. In order to make 
it k± dependent, an exponential term is added to get <f>^f Gen \ The third candidate DA model 
4>^ LS ^ was proposed in [ljj, which is a Gaussian type model. The fourth one [it] ], which is an 
exponential model, and the fifth model []j|], which is obtained by solving the equations of motion 
without three-parton contributions, were first proposed for B meson. Here we use heavy quark 
symmetry and modify the parameters to make them D meson DAs. The sixth DA was proposed in 
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FIG. 2: color suppressed diagrams in pQCD approach for B — > DP decays 
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FIG. 3: annihilation type diagrams in pQCD approach for B — > DP decays 
18|, which is derived from the BHL prescription [l^ j. with m<; = 0.35 GeV, m c = 1.3 GeV. In the 



above candidate DAs, only the second model has two parameters, and only (j)^ en ^ and (j)^ uan9 ^ 
are b independent. In the next section we will try to fit out the best D meson wave function 
parameters with the experimental results. As for D* meson, we just assume that 4>j^ = (fj^ = <p£) 
according to heavy quark symmetry. 



III. CALCULATION OF DECAY AMPLITUDES IN PQCD APPROACH 
A. Amplitudes for £>( s ) — ► D^P decays 

There are three types of diagrams that may contribute to the B — > D^*'M decays: color 
allowed diagrams (we mark this kind of contribution with the subscript ext) shown in Fig. 
color suppressed diagrams (marked with int) shown in Fig. [21 and the annihilation type diagrams 
(marked with exc) shown in Fig. [3j And each type of diagrams contains two categories: the one 
with one meson can be factorized out (denoted as £) and the one with no meson can be factorized 
out (denoted as M). 

The first two diagrams in Fig. Q3 [2] and [3] involve only two meson wave functions, whose results 



9 



are as following: 



/•l /-l/A 

^cxt = 8irC F f P / dx 1 dx 2 / bidb 1 b 2 db 2 (f)B(xi,b 1 )(f) D (x 2 ) 
Jo Jo 



x 



£ e (4 1} )/i(xi, x 2) 61, 6 2 )^(x 2 )(l + x 2 + r) + rE e (t^)h(x 2 , Xl ,b 2 , &i)St(xi)l , (31) 



6nt = 87rC F f D / (ixi(ix3 / bidb 1 b 3 db 3 (p B (xi,b 1 ) 
Jo Jo 

x{[(2 - x 3 )^p(x 3 ) - r (l - 2x 3 )(^(x 3 ) - 0?(x 3 ))] 
x^^ 1 ))^!, (1 - x 3 )(l - r 2 ), 61, &3)5 t (x 3 ) 
+2r ^(x 3 )^(tf - x 3 ,xi(l - r 2 ),b 3 , &i)S t (xi)} , 

?cxc = SttCf/b / (ix 2 da; 3 / b 2 db 2 b 3 db 3 <j) D (x 2 ) 
Jo Jo 



(32) 



x 



^30p(^3)-E'a(ti 1) )^a(2;2,X 3 (l - r 2 ) , 6 2 , 63)^ (x 3 ) 

+z 2 </>p(x 3 ).E a (4 2) )M ;z; 3, ^2(1 - ?" 2 ), h, b 2 )S t {x 2 ) 



(33) 



with the mass ratio tq = mo/ms- fp, fp and /d are the decay constants of the light pseudoscalar 
meson, B meson and D meson respectively. And the factor evolving with the scale t are given by: 

E e (t) = a s (t)ai{t)exp[-S B (t)-S D (t)], 
Ei(t) = a s {t)a 2 (t)exp[-S B (t) - S P {t)} , 

E a (t) = a s (t)a 2 (t)exp[-S D (t) - S P (t)] . (34) 

We adopt the expression of Sudakov factor for D meson as suggested in ref. Q], which is listed in 
appendix together with the expressions for Sp(t), Sp(t). 

The functions h's in the hard part of factorization formulae, derived from the factorizable 
diagrams, are given by 

h(x 1 ,x 2 ,b 1 ,b 2 ) 



h a (x 2 ,x 3 ,b 2 ,b 3 ) 



Kq (yJxix 2 m B bi) 
x [9(bi - b 2 )K {y/x^m B bi) I {y/x^m B b 2 ) 
+0(b 2 - h)K {y/x^niBbi) I (Vx^mph)] 



(35) 



i-j (v 'x 2 x 3 m B b 2 ) 

0(b 2 - b 3 )H^ {^/xim B b 2 ) Jo {\fxim B b 3 ) 
+9(b 3 - b 2 )H^ (y/ximsh) J {y/ximpb^ 



(36) 
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where H^ l \z) = Jq(z) + iYq(z). The hard scales t are determined by 

t*P = max(^/x^m B , I/&2) , t)p = max( v / i7m s , l/&i, I/62) > 



max 



(v/(l-x 3 )(l-r 2 )m B , 1/61,1/63) , *) 



(2) 



max 



(v/^i(l-r 2 )m B , I/61, 1/63) 



= m£ 



(Vx 3 (l-r 2 )m B , I/62, 1/63) , 4 2) = max( y^l - r 2 )m B , l/6 2 , 1/63) 



(37) 



The formulae for the last two diagrams in Fig. [H Fig. [2] and Fig. [3] contain the kinematics 
variables of the three mesons. Their expressions are: 

16iry / 2N c C F [ [dx] [ bidb 1 b 3 db 3 4>B{x 1 ,b 1 )(j) D (x 2 )(j)p(x 3 ) 
Jo Jo 

x x 3 E b (t^)h^(x i ,b l )-(l-x 3 + x 2 )E b (tP)h^(x i ,b i )\ , (38) 

1^^/W C C F [ [dx] [ 1 b 1 db 1 b 2 db 2 <t>B{xiM)<i>D{x2) 
Jo Jo 

((x 3 - 1 - x 2 )M*3) + r (l - x 3 )(<%(x 3 ) - f P {x 3 ))) E d {t^)h^\xiM) 



M 



+ [(1 -x 2 )(j ) p{x 3 ) +r (x 3 - l)(4> P P (x 3 ) + <$(s 3 ))] E d {tf)hf( Xi M 

\§-k^/W c C f I [dx] [ ' 6id&i6 2 d&20 B (xi,6 1 )<Mx 2 ) 
Jo Jo 

x 3 Mx3)E f (tf ) )hf\x i ,b l ) ~ x 2 Mx3)Ef(t { p)hf\xiA)] , 



(39) 



(40) 



with [dx] = dx\dx 2 dx 3 . The expressions for the evolution factors are 

E b (t) = a s (t)^-exp[-S(t)\ b2=bl ], 

JSc 

E d (t) = a s (t)^-e W [S(t)[ b3=bl ], 
E f (t) = a s ( t )^-ex P [-S(t)\ b3=b2 ], 



(41) 



with the Sudakov exponent S = Sp + Sd + Sp. 
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The functions j = 1 and 2, in these amplitudes are 



,00 



,00 



[6(h - b 3 )K (Bm B h) I (Bm B b 3 ) + 6{b 3 - h)K (Bm B b 3 ) I (Bm B h)} 
K (Bjm B b 3 ) 



for B] > 



fH^\^\B]\m B b 3 ) iovB]<0 
[0(&i - b 2 )K (Dmsh) I (Dm B b 2 ) + 9(b 2 - h)K (Dm B b 2 ) I (Dmsh) 
Ko(Djm B b 2 ) 



for D] > 



x < 



(42) 



fH^(^\D]\m B b 2 ) for D] <0 
if = i*- [fl(6i - b 2 )H^ (Fmsh) Jo (Fm B b 2 ) + 6(b 2 - h)H^ (Fm B b 2 ) J (Fm B h) 

Ko(F jmB h) for if > 
^ fH( 1 \ s f\Ff\m B b 1 ) for if <0 

with the variables 



(43) 



(44) 



S 1 



£ 2 2 



0? 



Dl 



F 



F- 



F. 



X\X 2 , 

x\x 2 - x 2 x 3 (l - r 2 ) , 

x\x 2 - x 2 (l - x 3 )(l - r 2 ) , 

xi(l - x 3 )(l - r 2 ) , 

(xi - x 2 )(l - x 3 )(l - r 2 ) , 

(zi + x 2 )r 2 - (1 - xi - x 2 )(l - x 3 )(l - r 2 ) , 

x 2 x 3 (l - r 2 ) , 

x 2 (xi - x 3 (l - r 2 )) , 

l-(l-x 2 )(l-x 1 -x 3 (l-r 2 )) . 



The scales £v') are given by 



4 j) = max(Sm B ,J|S?|mB, 1/63), 



4 j) = max(Dm B , J|D?|roB,l/&i,V&2) , 



4 j) = max(Fm B , A /|F?|m B , 1/62) • 



(45) 



(46) 
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The decay amplitudes of each BfA — ► D( S )P channels are then 

A (B~ -> D\-) = ^V cb V: d (Ut + M ext + Cint + M mt ) , (47) 

A (B- -> Z) ^-) = ^| (£ext + M earf + Cint + M int ) , (48) 

4 (50 _ £, +7r -) = + M ext + cj exc + M exc ) , (49) 

A (5° -> D+JT-) = (&*t + Mext) , (50) 

A (5° -> = (£exc + Mexc) , (51) 

A (B° - r>V) = ^|y cfe ^-L (_(^ nt + X . nt ) + + ^ exc )) , (52) 

A (B° - = ^=K*C + Mi^) , (53) 

4 (5O ^ D ^ n ) = ^L Vcb V* d ^= (Z in t + M M + £exc + M exc ) , (54) 

A {B° s - D+Tr-) = ^=K**& + AW) , (55) 

A (4° - ^tt-) = ( ^j=V d bVud (Ut + M ext ) , (56) 

A (4° -> ^ S + ^~) = ^| Vrf,C (£ext + M exf + Cexc + Mexc) , (57) 

A {B° s - L>V) = ^V cb V: s -L {Uc + Mexc) , (58) 

A {B° s - D i^ ) = ^VcbV* d (Ut + Mint) , (59) 

A (4° - £>V) = fexc + M exc ) , (60) 

A (4° - D% s ) = ^V cb V: s & nt + Mim) • (61) 

It should be noticed that, in (|54l) . (j60j) and (f6Tj) . the decay amplitudes are for the mixing basis of 

77 and rf . For the physical state r\ and rf , the decay amplitudes are 

A (5° D ??) = ,4 (B° -> D°rin) cos , (62) 

A (B° -» £> r/) = ^ (5° -» lA ?n ) sin0 , (63) 

A (B° -> D°r/) = A (5° D V) cos - j4 (5° D°r/ S ) sin , (64) 

A (B° s -> D°r/ ) = ^ (£° D V) sin + ^ (#° £> V) cos . (65) 
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B. Amplitudes for £?( s ) — ► D^V and B^ — ► DJ^P decays 

For the channels of i?( s ) — ► -D( S )V and D^P, no transverse polarization of the vector mesons 
will contribute. In the leading powder contribution, the formulae of B — > DV and B — > D*P are 
the same as that of Br s \ — > D( S )P decays, except some substitutions. 

For Br s \ — ► Df s \V, the following substitutions should be done for the formula & and A4f. 

(f>p -> 4>v , tf'p^ ~<t>v > 4>p -> ~<t>v ' rn ^m v , /p -> /y . (66) 

0y, 0y and </>y are the light cone distribution amplitudes of vector mesons, which we defined 
before. And my, fy are the mass and the decay constant of the vector meson related. 
For B( s ) — > D^P, the substitutions in the formula £j and A4j should be done as: 

m D -> m D » , /d -> f D * , c/>d(x 2 ) ->• 0d*(^2) • (67) 

Making the following substitutions in equations (|47p ~(|65j). we can get the final decay amplitude 
for each B — > D*P decays: 

D + _^ D *+ ; _^ D *o ^ D + _^ D *+ _ ^ 68 ^ 

And the formulae for B — > can be obtained through the substitutions 

7r -> p , if -> K* , ?? n -> uj , r? s ->• ^ (69) 

in equations (|l"T|) - (f61"j) . 

C. Amplitudes for B^ — ► £*( s )^ decays 

For S — > Df^V decays, both longitudinal and transverse polarization can contribute. For 
the longitudinal polarization, the amplitudes can be obtained by carrying out the substitutions 
referred in (|66p and (|67p , when only the leading power contribution is taken into consideration. The 
transverse polarized contribution is suppressed by r or ry, ry = my/ms- Although the transverse 
polarization will not give the leading power contribution, we still list the analytic formulae for 
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tT 

5CXC 



Scxt 



Sint 



■1/A 



transverse polarization £g xt , £ int and J 

87rCFm%fv / dx\dx 2 \ 6i(i&i6 2 (i& 2 <^ B (xi, 6i)<^(x 2 ) 
Jo Jo 

E e (tP)h(x 1 ,x 2 ,b 1 ,b 2 )S t (x 2 ) (-e n ^o^ T - ie *T . e *T {1 + 2r) 



+ J E e (4 2) )/i(x 2 ,xi,6 2 ,6i)5 t (x 1 )r(r + 1) 



. *T . *T 
%e D e V 



,1 ,1/A 

8itCFm B fD* \ dx\dx 3 / bidbib 3 db 3 <pB{xi,bi 



X /' 



^ ^ (-^(xa) - M(x 3 - 3)#(x 3 ) + (x 3 - l)^(x 3 ))) 
-ie*T . e f (^(jcg) _ ^((^ _ l)^(x 3 ) - (a* - 3)^(x 3 ))))] 
x^(tf ^(xi, x 3 (l - r 2 ), 61, 6 3 )S t (x 3 ) 



+7V 



e S«v (0« ( X3 ) - f v ( X3 )) + i e f • e* v T (r v (x 3 ) - <j> v v (x 3 )) 



xEi(tf >)% 3 , xi(l - r 2 ), 63, 6i)S t (xi) 



(70) 



(71) 



8-KC F m B $ B I dx 2 dx 3 I b 2 db 2 b 3 db 3 (^(x 2 ) 
Jo Jo 

x J B a (4 1 ))/ la (x 2 ,x 3 (l - r 2 ),6 2 ,6 3 )5 t (x 3 ) 

x ( e nne^^ [r 2^ (2 , 3) _ + + rry(1 _ x 3 )^(x 3 )] 

+ie* D T ■ ef [-r 2 ^(x 3 ) + rr y (x 3 - l)^(x 3 ) + rr v (x 3 + l)^(x 3 )]) 
+ E a (tW)h a {x 3 , x 2 {l - r 2 ),b 3 ,b 2 )S t {x 2 )rr v 
x (e nSe ^ T ((l + x 2 )<^(x 3 ) + (1 - x 3 )^(x 3 )) 

-ief ■ 4 ((1 - x 2 )^(x 3 ) + (1 + x 2 )^r(x 3 )))] • (72) 

The evolution factors in the amplitudes are the same as those in (J35]) and (|36|) after substitut- 
ing Sy(t) for Sp(t). For the nonfactorizable amplitudes, the factorization formulas involve the 



15 



kinematic variables of all the three mesons. Their expressions are 



M 



cxt 



1/A 



16iry/2N c C F m% / [dx] / b 1 dbib 3 db 3 4> B {x 1 ,b 1 )(p] D (x2) 



E h {t^)hf\ Xi M 



x (e nne ° *v x 3 (^(x 3 ) - <fy(x 3 )) + ie£ ■ x 3 (#(x 3 ) - (/> v v (x 3 )) 
+^ 6 (4 2) )^ 2) (x il 6i) 

x {e nn *° e v((l - x 3 )(l - 2r)<^(x 3 ) + (x 3 - l)^(x 3 )) 
-ze* D T • ef ((x 3 - l)^(x 3 ) + (1 - 2r)(l - x 3 )^(x 3 ))} 

A*£ t = lforv^Cf-mB / [dx] / 6id&iM&2<Mxi,&i)^(x 2 ) 

jo Jo 



xr 



E d {tf)hf{xiM 



x e 



4 (( X2 - l)<^(x 3 ) + M<^(x 3 ) + <^(x 3 ))) 



-if 



f ((1 - x 2 )^(x 3 ) + ry(^(x 3 ) + <^(x 3 )))) 



At 



+^# ) )/# ) (*i,&0 

X (e^^r^X, - 1)^(X 3 ) - X 2 ^(X 3 )) 

-ie*S ■ e* v T (x 2 (P t v (x 3 ) + 2r y (x 3 - l)^(x 3 )))] 

1/A 



(73) 



(74) 



= l6ir^2N c CFm% / [dx] / b 1 dbib 2 db24>B{xi, 6i)^(x 2 ) 
Jo Jo 

X {^(^^A^v^s) 

x ( e ^^(x 2 r 2 - 4x 3 ) - fc^ • ef(x 2 r 2 + ^x 3 )) 
+S / (4 2) )4 2) (x,,6 i ) 

x " e «-B T ^ T (_2rr y ^(x 3 ) - r 2 (x 2 - l)^(x 3 ) + 4(x 3 - l)^(x 3 )) 

+^*D T • (((x 2 - l)r 2 + r 2 v (x 2 - l))^(x 3 ) + 2rry^(x 3 ))] | . (75) 

The fo's and /&^°) functions in the amplitudes here are the same as defined in (|35l36l42M44p , 

If carrying out the substitutions in (f68j) and ([69]) together, we can get the final decay amplitude 
of each polarization for B — ► D*V decays. 
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TABLE II: The smallest \ 2 for each kind of the D meson DAs 



$(Gen) $(MGen) §(KLS) §(GN) §{KKQT) §(Huang) 

34.1 33.6 156.9 112.3 45.0 41.9 



IV. NUMERICAL CALCULATIONS AND DISCUSSIONS 

(*) 

In this section the numerical results of our calculation will be given. The parameters of 
meson we use are 

m D - = 1.869GeV, m n - = 1.968GeV, 
m D *- = 2.010GeV, m n *- = 2.112GeV, 
f D - = 223MeV, f D - = 274MeV. (76) 

A lot of study has been made on the decay constants of Dr s \ mesons. Here we take the values from 
ref. [2^]. Since there are no experimental results of the decay constants of D* s -^ mesons, we use the 
relations between fo and fo* derived from HQET: 



fn*- = J^fn-, f Dr = J^-f D -, (77) 



mo*- \ mo* 



which is different from ref. 



2l|. 



With the D meson wave functions we get, we can calculate the amplitudes easily. And the 
decay width can be got by 

r = ^G 2 F \V cb \ 2 \V ud \ 2 m 7 B (1 - r 2 )\A\ 2 , (78) 
with A the decay amplitude defined in eqs. (|47H6"T1) The branch ratio is 

Br = Th/T Bw , (79) 

with tb (s) as the life time of meson. We take r B - = 1.674 x 10~ 12 s, t b q = 1.542 x 10 _12 s, 
t b o = 1.466 x l(T 12 s, and G F = 1.16639 x 10" 5 . 

A. results of fitting 

Since the B — > DP decay channels have been measured in high precision, we use these ex- 
perimental results to fit out the parameters of the candidate D meson DAs. Here we don't use 
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FIG. 4: The D meson distribution amplitude <f)p Gen (0.8, 0) ( the blue, solid line) and D s meson distribution 
amplitude </>£f Ge ™ (0.5,0) (the red, dotted line) 



the experimental results containing r] or rj' in the final states because there are uncertainties from 
the mixing. The six decay channels we used to fit out the D meson wave function parameters 
are B~ -> D°ir~, B~ -> D°K~,B° D+ir~, B° -> D + K~, B° -> L>°vr , 5° -> L>°if°. The 
experimental results of these channels are listed in table IIIH which are from ref. 
we used for fitting is 



221. The formula 



X 



£ 



07 



(80) 



The i means the summation over the six decay channels. Br? x (Brj h ) is the experimental (theoret- 
ical) value of branch ratio, and Oj is the uncertainty of the experimental value. In table [TT] we list 
the smallest x 2 we get f° r an the D meson DAs. Easy to see, except 

$(kls) and $(GAT) ) all the 

other DAs have a small Xmin- &( MGen ) is the best one, with its parameters fixed as Cr> = 0.8, 
ijj = 0.1. We will use this D meson wave function for the following numerical calculations of all the 
decay channels. For the D s meson, we use Co = 0.5, u = 0.2, with a little SU(3) breaking effect. 
In this case, we can see from Fig. U that the s quark in D s meson has a little larger momentum 
fraction than the d quark in the D meson, which characterize the little larger mass of s quark. 
Because the mass difference between the vector meson , and pseudoscalar meson D( s ) is 



small, so we adopt the same DA for them. 
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TABLE III: Branching ratios of Br s ^ — ► DP decays calculated in pQCD approach with experimental data 
(in units of 10~ 4 ) 







Experimental results 


Our results 


B~ - 




47.5 ± 1.9 


ci -,+29.5+4.3+1.5 

— Z\J. { — ( .O— l.O 


B~ - 


-> D°K- 


3.83 ±0.45 


a nn+2. 35+0. 63+0. 12 
^^ — 1.64—0.93 — 0.12 


B°- 


> L»+7T- 


26.5 ± 1.5 


9n Q +17. 8+5.5+0.8 

— 11. ( — ( .0 — U.o 


B°- 


> £>+A'- 


2.04 ±0.57 


9 ^ o + l. 56+0. 63+0. 07 

— 1.01 — 0. ( 1 — U.O ( 


B°- 


> Z?°7r° 


2.61 ±0.25 


1 QC + 0. 67+0. 51 + 0.06 
— U.OD— U.OO — U.UO 




> D°K° 


0.523 ±0.066 


n oo +0.08+0.06+0.01 

— U.Ui — U . UD — U . U 1 


B°- 




2.02 ±0.21 


9 /j^+0. 97+0. 30+0. 07 
— u. / / — u.ou — U.Ui 


B°- 


> L>y 


1.26 ±0.21 


1 fic+0. 65+0. 20+0. 05 
— U.o^ — U.zU- U.Oo 


B° - 




269 ± 054 


u -^ o -0. 12-0. 11-0. 01 


B.°- 


-> D+7T" 




9 O9+0.82+0.32+0.07 -, n -2 
z - oz -0. 61-0. 39-0. 07 x lu 


B°- 


Z? rr° 




1 i c+0.36+0. 19+0.04 1n -2 
x -- Ld -0.29-0.20-0.04 x lu 


B°- 


Z? A° 




o QC + 1. 49+0. 88+0. 11 
°- yo -l. 25-0. 99-0. 12 


B s °- 






n i 1+0.05+0.03+0.00 
u - i *-0. 05-0. 03-0. 00 


B s °- 






n oo + 0. 11 + 0.03 + 0.01 
u - oz -0. 10-0.04-0.01 


B s °- 




38±3±13 


91 o + H.4+6.9+0.6 
zi '°-8.1-6.8-0.6 


B s °- 






1 71 +0.92+0.58+0.05 
L - 1 i -0. 65-0. 55-0. 05 



B. Results for all the related channels and discussions 

Our numerical results are listed in table IIIIl table IIVI table [V] and table IVII The first error in 
these entries are caused by the hadronic parameters in B/ s \ meson wave function (the decay constant 
and the shape parameter). We take f B = 0.19 ± 0.025, co b = 0.40 ± 0.05 and f B , = 0.24 ± 0.03, 
uj? s = 0.50 ± 0.05. The second error arise from the higher order perturbative QCD corrections: 
the choice of the hard scales, defined in (137p and (|46p . which vary from 0.75t to 1.25t, and the 
uncertainty of ^\qqd = 0.25 ± 0.05. The third error is from the uncertainties of the CKM matrix 
elements. In our calculation, we use 

V cb = (41.6li°f) x 10- 3 , V ud = 0.973851°;°^ , V us = 0.22715^°;°oioo ■ (8 1 ) 

Among them, the hadronic inputs always gives rise to the largest uncertainty, and the CKM matrix 
elements contribute little. 

The first six channels in table IIIIl are input values of the x 2 fit program. Although we get a 
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TABLE IV: Branching ratios of B^ — > DV decays calculated in pQCD approach with experimental data 
(in units of 10~ 4 ) 







Experimental results 


Our results 


B~ - 


■+ D°p- 


134 ± 18 


i i o+68. 1+18.8+3.3 
^—45.9 — 26.1 — 3.3 


B~ - 


-> D°K*- 


5.29 ±0.45 


p jq+3. 86+0. 12+0. 20 
— z.oo — l.oo— U.zU 


B°- 


> D+p- 


75 ± 12 


an o+45. 0+13. 1+2.0 
— zy . i — i (.© — z.u 


B°- 


> D+K*- 


4.60 ±0.78 


a H7+ 2 - 61+0. 94+0. 12 
— i.oy — i.ii — u.iz 


B" - 


>D°p° 

r 


2.91 ±0.40 


i 7Q+0. 65+0. 46+0. 05 
— U.o2— U.oy— U.05 


B° - 


>D°LJ 


2.60 ± 0.29 


a r + l. 58+0. 85+0. 12 

— 1 . Zi — U. io — U. lz 


B°- 


> D°K*° 


0.423 ±0.064 


n o/?+0. 10+0.08+0.01 

— U . Uo — U.U 1 — U . U 1 


B" - 


> D + K*~ 


< 8 


i qa+0. 66+0. 41+0. 06 


B s °- 






n i i +0.03+0.02+0.00 
u ' 11 -0. 03-0. 02-0. 00 


B s °- 


D°p° 




(r , 9 +1.57+0.96+0.16\ 1n -2 
V -* -1. 39-1. 11-0. 16/ x iu 


B s °- 


-> D°K*° 




A co+1.70+1.15+0.13 
*- dz -1.43-1.30-0.13 


B s °- 


■4 D°LJ 




(a 71 +1.35+0. 88+0. 14\ v 1n -2 
V±-' 1 -1.18-0.92-0.14j x 1U 


B s °- 






n qn+OU+0.07+0.01 
u - ou -0. 10-0.08-0.01 


B s °- 






ci n +26. 9+16. 6+1.48 
— 19.2 — 16.2 — 1.48 


B s °- 






q no+l. 62+0. 88+0. 10 
°- uz -l. 16-0.91-0. 10 



reasonable x 2 m the fit, the branching ratio of B° — > D°K° is about only half of the experimental 
value. Comparing with the color suppressed diagrams, the annihilation diagrams contributes little 
for B° — > D°K° and B° — > D°7r°. So we can take them for a compare. They have different CKM 
elements {V c bV* s for the former and V c t>V* d for the latter). Taking the factor -4^ in flavor wave 
function of ir meson into account, the Br (B° — ► D°K°) is roughly one tenth of Br(B° — ► D°7r°). 
So this value of Br(S° — ► D°K°) is theoretically reasonable. The similar argument is valid for 
Br(S° -» 

Although we use only six B — > DP channels to fix the D meson wave function, the results 
of most other channels especially those B — > DV and D*P channels agree very well with the 
current experimental measurements. Easy to find that, Br(i? — ► DoS) is twice larger than Br(i? — > 
Dp), while their experimental results are near to each other. Both these two channels receive 
contributions from the color suppressed diagrams and annihilation diagrams. They are at the 
same order magnitude for these two processes. For color suppressed diagrams, the dd of the flavor 
part contributes. While the uu part contributes to the annihilation diagrams. Amplitudes of these 
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TABLE V: Branching ratios of J3( s ) — > D*P decays calculated in pQCD approach with experimental data 
(in units of 10~ 4 ) 



Experimental results Our results 



B~ - 




52.8 ±2.8 


r n 4 +29.4+4.4+1.5 
ou -^-20.4-7.3-1.5 


B~ - 


-> D*°K~ 


3.6 ± 1.0 


o qo+2. 33+0. 62+0. 12 
°- y °-l. 62-0. 92-0. 12 


B°- 


* D*+n- 


26.2 ± 1.3 


OR n+ 17 3 + 5 - 34 + - 75 
zo - u -ll. 4-7.09-0. 76 


B°- 


* D*+K~ 


2.04 ±0.47 


n 07+1. 52+0.62+0.07 
z -°' -0.99-0.69-0.07 


B a - 


> D*°7T 


1.71 ±0.28 


9 1 1 +0.76+0.52+0.06 
<fi.J-J._0. 69-0.62-0. 06 


B°~ 


* L>*°r/ 


1.80 ±0.31 


9 Kn +1. 01+0. 29+0.08 
z - DU -0. 83-0. 33-0. 08 


B°- 


> £>* ?/ 


1.21 ±0.40 


1 7/1+0. 68+0. 19+0. 05 
1 - '^-0.55-0.22-0.05 


B°- 


> D*°K° 


0.36 ±0.12 


n 9/1+0.09+0.06+0.01 
u - z ^-0. 07-0. 06-0. 01 


B a - 


* X^+A'" 


0.200 ±0.064 


n 47+O.I6+O.I2+O.OI 
u -^' -0.12-0.11-0.01 


B s °- 


-> A>*+7r- 




jr) 97 +0.81+0. 31+0. 01 \ 1n -2 
l z - z ' -0.60-0.38-0.01/ x iu 


Bs°- 


-> D*°7r° 




/ 1 1 o+0.35+0.18+0.03\ v 1n -2 
V i - ±o -0. 28-0. 19-0. 03/ A iu 


B s °- 


-> D*°K° 




/, 97+I.54+O.88+O.I2 
*- z '-1.35-1.05-0.12 


B s °- 


-> £>*°r/ 




n 1 c+0. 06+0. 03+0. 00 
u - io -0. 05-0. 03-0. 00 


B s °- 


£>*V 




n oo+0. 12+0. 03+0. 01 
u -°°-0. 10-0.03-0.01 


B s °- 


-> A>:+7r- 




94 9 + 11.2+7.8+0.7 
&±.&_rj .2-7.7-0.7 


B s °- 


-> A»*+A- 




1 GC+0.90+0.56+0.05 
J" od -0.63-0.53-0.05 



two kinds of diagrams have the same sign in the B — ► Du> decay but different sign in the B — > Dp 
decay due to isospin. Similar situation exists for Br(S — ► D*uj) and Br(l? — > D*p). 

The B° — ► DfK~ decay is a kind of pure annihilation type decays dominant by W exchange 
diagram. Our result is larger than the experiments and also larger than the previous pQCD 
calculations [23|] due to the change of D meson wave functions. The annihilation type diagrams 



are power suppressed in pQCD approach, which is more sensitive to the hadronic wave functions. 

For the decays B — > in table IVIl we also estimate the ratios of transverse polarized 

contribution Rt = \ At\ 2 / (|^4t| 2 ± |^4l| 2 )- We should mention that, these results are just indicative, 
because transverse polarizations are power suppressed by ry or r to make it more sensitive to 
small parameters and higher order contributions than the longitudinal contribution. Although the 
transverse polarization is suppressed in B — > D*V decays, in some channels, such as B — ► D*°p° 
and B -> D*°K*°, etc, it has 40% contributions. The reason is that the dominant contribution in 
these channels is from Mint in eq. H39|) . which is X3 suppressed, while the transverse contribution 
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TABLE VI: Predicted branching ratios of Br s -) — » D*V decays with experimental data (in units of 10 4 ) 
together with the percentage of transverse polarizations Rt 







Experimental BRs 


BRs in pQCD 


Rt 


B- - 


■4 D*°p~ 




i i 7+71.0+19.9+3.4 
11 ' -48.0-27.0-3.4 


0.04 


B~ - 


-> D*°K*- 


8.3 ±1.5 0.14 


a O9+4.14+1.22+0.21 
u ' oz -2. 80-1. 65-0. 21 


0.06 


B°~ 


> D*+p- 




70 ,(+52.6+15.9+2.3 
' y '^-34.5-21.2-2.3 


0.15 


B°- 


> D*+K*- 


3.20 ±0.67 


/I oo+3. 18+1. 16+0. 15 
*-oo_2. 08-1. 41-0. 15 


0.19 


B°- 


> D*°p° 


3.73 ±0.99 


q RQ+1- 52+0. 57+0. 11 
O.U»_ 1 23-0 .63-0 .11 


0.47 


B°- 


>D*°to 


2.68 ±0.50 


r oq+2. 14+0.84+0. 17 
84-1. 10-0. 17 


0.24 


B°- 


> D*® K*® 


< 0.69 


n m +0.20+0.08+0.02 
u - ol -0. 17-0. 10-0. 02 


0.44 


B°- 


> D*+K*- 




i Q7+0. 58+0. 37+0. 06 
1 - y 1 -0.49-0.43-0.06 


0.02 


Bs°- 


-> D*+p- 




n i i +0.03+0.02+0.00 
u - ii -0. 03-0. 02-0. 00 


0.01 


B°~ 


-> D*°p° 




(r nc+1.44+0.96+0.16\ i n -2 
l°- 00 -l. 39-1. 05-0. 16 J x iu 


0.01 


B s °- 






o ao+3. 30+1. 68+0. 24 
o-4d_2. 67-1. 83-0. 25 


0.43 


B s °- 


-> D*°uj 




(A fi n+l-50+0.89+0.14\ lf) -2 
l^- DU -l. 11-0.84-0. 14j x iu 


0.01 


B s °- 






n ri +0.19+0.09+0.02 
u - oi -0. 17-0. 11-0.02 


0.37 


B s °- 






c R q+30. 4+19. 1+1.65 
<-""•y_21.6-18.3-l.66 


0.13 


B s °- 


-> D*+K*- 




o 4 7+1.96+1.07+0.11 
-1.35-1.06-0.11 


0.17 



in eq, (|74j) is only r suppressed. They are comparable numerically to make a large contribution for 
transverse polarizations in these color suppressed channels. This mechanism is different from those 
charmless B decays where the dominant transverse polarizations are from the space like penguin 
(penguin annihilation) contributions 



271 ]. Here the annihilation type contributions are mainly from 
W exchange diagrams contributing little to transverse polarizations. 

For those previous calculated channels in pQCD approach 0,2], our results are slightly different 
due to parameter changes. Most of the B (B^) decay channels are measured by the two B factories, 
which are consistent with our calculations. For the B s decays, only one channel is measured. Our 
predictions will soon be tested by the LHCb experiments. 

For comparison with other methods, we also give the form factors at the maximal recoil 

cB^D _ n c 9 +0.15+0.05 tB s ^D s _ n ,,,+0.11+0.07 ( on\ 
^ + — U.O__ 012 _ 07 , ^ + — U. _O_ 09 _ 08 . (0_j 

These are comparable with other methods 2a] . 
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If applying the naive factorization approach, we can get 

A{ B° - D+tt") = i^L Vcb V: d {M 2 B - M 2 D )UF B - D {Ml) ai {D^), (83) 
V2A(B° - DV) = -i^ Vcb V: d (M 2 B - Ml)f D F B ^(Ml)a 2 (DK). (84) 

Substituting our results for A(B° -> D + ir~) and ,4(5° -> L>°vr ) in eq.® and ([84]) . we can 
extract the BSW parameters a\ and 02 from our pQCD approach 

\a 2 /ai\ = 0.49 , Arg(a 2 /a 1 ) = -37.6°. (85) 

If the annihilation diagrams' contribution is excluded, the results are 

\a 2 /ai\ = 0.54 , Arg(a 2 / ai ) = -59.0°. (86) 

Indeed, the large | a 2 /ai | implies that the color suppressed decays not very suppressed as previous 
expected Q]. The relative strong phase between the two contributions is not small ^as naive ex- 



pectations. These results are consistent with recent direct studies from experiments 
difference is that our results come from direct dynamical calculation not from fit. 



2J. But the 



V. CONCLUSION 



In this paper, we calculate the branch ratios of —>■ D^P, D*^P, D^V and D*^V channels, 
with the D meson wave function obtained through fitting. We also calculate the ratios of transverse 
polarized contributions in B —> D*V decays. Most of the results agree well with the experiments. 
It seems that there's a disagreement with the experimental data in the relative size of branching 
ratios for B D^p and B -> £>(*) u. Some channels of the B —> D*V decays may receive large 
contribution from the transverse polarization. The results of B® — > D^P, D^V, D* s ^P, D* s ^V 
decays will be tested in the future experiments. 
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APPENDIX A: PQCD FUNCTIONS 



Jet function appears in (|3B - p3l . ([70J)-([72]) is 



(Al) 



the value of c in the above equation is 0.5 in this paper. And the Sj(xi)(j = B,C,P or V) functions 
in Sudakov form factors in (jM)) and ([^T|) are 

"* dp, 



S B (t) = s [ xx^,h ) +2 / ^ q (a s (fj,)), 

S c (t) = s[x 2 —^,b^\ +2 I ^7 ? (o; s (/2)), 
V2 / J\/b 2 V 



(A2) 
(A3) 



Sv(t) = S P {t) = s^x^,b^+s(^l-x^,b^+2j t ^ ^ 7<z K(/2)), (A4) 
with the quark anomalous dimension y^ = —a s /ir. The explicit form for the function s(Q, b) is: 



AW 
2/3i 



A( 2 ) A« /e 2 ^" 1 
In 



" 4/3? 



ln(2g) + 1 ln(26) + 1 



+ 



8/3? 



4/3? 4/?! 
In 2 (2g) - In 2 (26) 



In I ± 

6 



where the variables are defined by 

g = ln[Q/(^A)], b = ln[l/(6A)], 
and the coefficients and $ are 

33 - 2n f 



(A5) 



(A6) 



01 



-, ft 



153 - 19n 



12 ' 24 
67 vr 2 10 



. ri A 4 67 10 8 .1 . 

A^ = -, A (2) = n# + -0i n -e 7B , 

3' 9 3 27 ' 3 M v 2 ; ' 



(A7) 



is the number of the quark flavors and je is the Euler constant. We will use the one-loop 
running coupling constant, i.e. we pick up the four terms in the first line of the expression for the 
function s(Q,b). 
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